Astronomy <fc Astrophysics manuscript no. h3507 


February 1, 2008 


(DOI: will be inserted by hand later) 





Efficient acceleration and radiation in Poynting flux powered 

GRB outflows 

G. Drenkhahn* and H. C. Spruit 

Max-Planck-Institut fur Astrophysik, Postfach 1317, 85741 Garching bei Miinchen, Germany 
Received 20 February 2002 / Accepted 4 June 2002 

Abstract. We investigate the effects of magnetic energy release by local magnetic dissipation processes in Poynting 
flux-powered GRBs. For typical GRB parameters (energy and baryon loading) the dissipation takes place mainly 
outside the photosphere, producing non-thermal radiation. This process converts the total burst energy into 
prompt radiation at an efficiency of 10-50%. At the same time the dissipation has the effect of accelerating the 
flow to a large Lorentz factor. For higher baryon loading, the dissipation takes place mostly inside the photosphere, 
the efficiency of conversion of magnetic energy into radiation is lower, and an X-ray flash results instead of a GRB. 
We demonstrate these effects with numerical one-dimensional steady relativistic MHD calculations. 
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1. Introduction 

High luminosity outflows from 7-ray bursts (GRBs) must 
have large Lorentz factor to overcome the compactness 
problem (e.g. Piran 1999). To fulfil this requirement the 
mass loading can only be small so that the total en- 
ergy density exceeds greatly the rest mass energy den- 
sity. Poynting flux can carry large energy amounts though 
vacuum which provides a mechanism to transport energy 
without the need of matter. The release of electromagnetic 
energy by the central engine of a GRB is part of many 
models. E.g. tori in merger scenarios may by highly mag- 
netised due to the field amplification by the differential ro- 
tation (Narayan et al. 1992; Thompson 1994; Meszaros & 
Rees 1997; Katz 1997). Alternative models involve highly 
magnetised millisecond pulsars (Usov 1992; Spruit 1999). 
In all cases the rotational energy of a compact object 
will be tapped and the rotating magnetic field produces a 
Poynting flux. 

While Poynting flux is thus a plausible way of power- 
ing a GRB, it is not a priori clear how this energy flux 
is converted into the observed 7-rays. To accelerate the 
matter to the observed high Lorentz factors, a part of the 
Poynting flux must be converted into kinetic energy. This 
energy later powers the afterglow when it is released in an 
external shock. Since the prompt emission in most cases 
accounts for the bulk of the observed radiation, a mech- 
anism is needed to efficiently convert a magnetic energy 
flux into non-thermal radiation. 
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For the acceleration of the flow one can think of mag- 
netocentrifugal effects. But trying to explaining the flow 
acceleration by stationary ideal MHD processes is prob- 
lematic. For the purely radial magnetically driven stellar 
wind (modelled originally by Weber & Davis 1967; Mestel 
1968) the radial gradient of the magnetic pressure and the 
inward pointing magnetic tension force act against each 
other. Especially in the case where the flow is relativis- 
tic from the beginning there is a balance between these 
forces so that such purely radial flows are not accelerat- 
ing. However, an acceleration takes place if the flow lines 
diverges faster with radius than in the radial case beyond 
the fast critical point. The tension and pressure gradient 
forces are out of balance and Poynting flux to kinetic en- 
ergy flux conversion occurs (Begelman & Li 1994; Daigne 
& Drcnkhahn 2002). A magnetic acceleration model which 
uses only ideal MHD thus has to provide just the right flow 
divergence for acceleration to take place. As we show in 
this paper, a better alternative is a flow in which part of 
the magnetic energy density is dissipated locally. The de- 
crease of magnetic energy density with distance in such 
a model causes effective acceleration (Lyubarsky & Kirk 
2001; Drenkhahn 2002), while at the same time providing 
an efficient energy source for the observed 7-rays. 

The currently most accepted model explaining the 
high energy prompt emission of GRBs is the internal 
shock model (Rees & Meszaros 1992, 1994; Sari & Piran 
1997). Variations of the central engine luminosity pro- 
duces flow shells with different Lorentz factors which col- 
lide. Through these collisions a part of the kinetic energy is 
transfered into prompt radiation. The energy conversion 
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is only efficient if the spread in Lorentz factors is large 
(Kumar 1999; Panaitescu et al. 1999; Beloborodov 2000; 
Kobayashi & Sari 2001). The observed ratio of afterglow 
and prompt emission indicates a high efficiency of energy 
conversion into prompt emission. While this observation 
does not rule out the internal shock model, it does put 
strong constraints on it. 

If one allows for non-ideal MHD processes magnetic 
energy can be transfered to the matter by dissipation 
through reconnection. For this to happen, there must be 
small scale variations in the magnetic field. The energy 
that is released by washing out these variations can be 
converted into radiation. We call this the 'free magnetic 
energy' in the flow. An example of such small scale varia- 
tions would be the 'striped' wind (Coroniti 1990) that re- 
sults from the rotation of an inclined magnetic dipole. The 
distance between neighbouring stripes of different field di- 
rection in this case is ttv/H,, where v is the outflow speed 
and VL is the dipole's angular frequency. This scale is of 
the order of the size of the central engine (assumed to 
be a relativistic object). In general, all non-axisymmetric 
components of the magnetic field of the central engine pro- 
duce such variations. If the magnetic configuration is pre- 
dominantly non-axisymmetric, almost all of the Poynting 
flux is in the form of a field that changes sign on such a 
small length scale. This is the model we will use for the 
quantitative calculations below. We note, however, that 
even an axisymmetric rotating field can in principle pro- 
duce small length scales. The outflow near the axis of an 
axisymmetric MHD flow is spiral-like. This configuration 
is kink unstable so that field components can reconnect 
across the rotation axis. For a discussion of this point see 
Spruit ct al. (2001), hereafter Paper I. As we showed in 
Drenkhahn (2002) (Paper II from here on), such kink- 
produced irregularities are somewhat less efficient at con- 
verting magnetic energy. Since perfect axisymmetry is a 
special case we regard the non-axisymmetric case in this 
study. 

Fast reconnection leads to a decay of the magnetic 
field. The flow accelerates since the field decay induces an 
additional outward gradient in the magnetic pressure. In 
Paper II we explored the dynamical effects of the Poynting 
flux dissipation in the outflow. We found there that for 
fiducial GRB parameters a large amount of the Poynting 
flux energy is converted to kinetic energy. Also, a great 
part of the dissipation happens in the optically thin re- 
gion so that a potentially large fraction could be converted 
into non-thermal, prompt radiation. This model offers an 
alternative to the internal shock model in explaining the 
prompt emission by local dissipation of free magnetic en- 
ergy. 

The results presented in Paper II were based on an 
analytic approximation for the flow. In the present paper, 
we relax this approximation, and analyse a Poynting flux 
powered wind numerically. The results confirm the main 
results from the analytical study, but in addition allow 
us to determine which fractions of the Poynting flux are 
converted into thermal and non-thermal radiation. In this 



way, we can also determine the conditions under which a 
true GRB, as opposed to an 'X-ray flash' (Heise et al. 2001; 
Heise & in 't Zand 2002) is produced by a magnetically 
powered outflow. 

2. The Model 

We consider a radial outflow of magnetised plasma with 
the magnetic field being aligned transversal to the flow 
direction. The field contains small scale variations in di- 
rection from which energy is released. We parameterise 
the variation by introducing a variation length scale A w 
27rc/f2 on which the direction of the field changes. Field 
variations in the outflow are naturally produced by any 
non-axisymmetric component of a rotating magnetic field. 

The rate of magnetic energy dissipation is governed by 
the reconnection rate between neighbouring regions of dif- 
ferent field line direction. For highly symmetric initial con- 
ditions, the initial reconnection process is sensitive to the 
microscopic diffusion rate, but this situation is rarely rele- 
vant in astrophysics. Instead, the field reconnects by 'rapid 
reconnection' processes, in which the reconnection speed 
depends only logarithmically on the microscopic transport 
coefficients (Petschek 1964; Parker 1979; Priest & Forbes 
2000). Since the reconnection rate is an important factor 
influencing the results, we keep track of its effect by a 
parameter study. For this purpose we write the reconnec- 
tion time scale r co across the variation length scale A co as 
T co = Aco/^r, where the reconnection speed v r is the veloc- 
ity at which field lines of different directions are brought 
together by the dynamics of the reconnection process. We 
regard this process in the comoving frame moving with the 
bulk large-scale flow. The speed v r is known to scale with 
the Alfven speed va, i.e. v r = sva where e is a numerical 
factor < 1. For rapid reconnection in 2 dimensions, for ex- 
ample, numerical results (e.g. Biskamp 1986) show that e 
can be of the order 0.1. Since we do not know the density 
of reconnection centres the overall rate of field dissipation 
is still unknown. By adjusting the parameter e towards 
lower values we compensate for this ignorance. Thus, we 
take e as a measure for both, the reconnection speed and 
the density of reconnection centres. 

As a second parameter of less importance, we introduce 
the fraction /i 2 of the magnetic energy density that cannot 
be dissipated by local reconnection. If the field is the re- 
sult of the winding-up of a completely non-axisymmetric 
field (i.e. with vanishing azimuthal average), the direc- 
tion of the field lines changes completely over one varia- 
tion length, and we have [i = 0. If, on the other hand, 
the axisymmetric component does not vanish, we have 
fi = B-ft/B, where and B are the amplitudes of the 
axisymmetric component and the total field. In most of 
our study, we will assume that the Poynting flux decays 
completely which correspond to fi = 0. 

The released magnetic energy is converted into ther- 
mal and kinetic energy. The direct conversion into kinetic 
energy is possible because the field dissipation induces an 
additional outward gradient in the magnetic pressure. So 
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even if there is no thermal pressure gradient (as in the 
cold approximation used in Paper II) part of the released 
energy accelerates the flow directly. 

Though small scale structures in the flow are an essen- 
tial part of the model, we will only work with quantities 
which are averaged over these small scales. One can con- 
sider the flow to be stationary on large length scales so 
that we can do a time independent calculation to obtain 
general results. The adjustable variable e accounts for the 
unknown processes on small scales and includes effects in- 
troduced by the averaging. 

The dissipated magnetic energy initially takes the form 
of internal (thermal) energy of the gas. If the cooling time 
is long compared with the expansion time scale of the flow, 
the energy is mostly converted to kinetic energy, through 
adiabatic expansion. This is the case in the optically thick 
part of a GRB flow, inside its photosphere. The small 
fraction of the thermal energy that remains when the flow 
passes through the photosphere then shows up as thermal 
radiation emitted at the photosphere. In the optically thin 
parts of the flow, the radiative cooling times are typically 
quite short compared with the expansion time scale as 
will be shown in Sect. 2.6. The medium stays cold, and 
all thermal energy gets quickly converted into radiation. 
We assume the radiation processes to be similar to those 
invoked in the internal shock model, so that the radiation 
is non-thermal. Thus, we can only estimate the spectrum 
of the (small) part that is emitted as thermal radiation at 
the photosphere. We are able to calculate the total amount 
of non-thermal radiation, however, since this only depends 
on the rate magnetic dissipation. 

We make a two-zone approach with respect to the op- 
tical thickness to simplify the treatment. The flow is op- 
tically thick up to the photosphere and matter and radi- 
ation are treated as one fluid in permanent thermal equi- 
librium. At the photosphere the internal energy carried 
by the radiation decouples from the matter and escapes 
as black-body radiation. From the photosphere on, part of 
the Poynting flux produces non-thermal radiation through 
the magnetic dissipation process, while the rest still accel- 
erates the flow. 

2.1. Evolution of the magnetic field 

We model the evolution of the magnetic field B in a 
Poynting flux dominated outflow by a dissipation time 
scale t which depends on the Alfven speed and a typical 
length scale of the field geometry considered. The moti- 
vation and detailed derivation for this approach is given 
in Paper II. Since the dissipation time scale r depends on 
the local Alfven speed in the flow it is a function of the 
proper mass density p, the proper internal energy e, the 
absolute value of the radial bulk 4- velocity u and the mag- 
netic field strength B: r = r(p, e,u, B). We consider all 
thermodynamic quantities p, e, . . . in the comoving frame 
while all other quantities B,u,t, . . . refer to the lab frame, 
the frame in which the central engine rests. We use the 



notation T = ^/\ + u 2 for the bulk Lorentz factor of the 
outflow and (3 — u/T for the bulk velocity in units of the 
speed of light c. 

The flow is assumed to be purely radial, and we are 
considering distances from the central engine that are suf- 
ficiently far from the Alfven radius. Thus any centrifugal 
acceleration of the flow has already taken place, and the 
dominant field component is £></, = B 3> B r , Bg . Without 
internal dissipation, the induction equation would thus 
yield d r (3rB = 0. The evolution equation for the mag- 
netic field including dissipation is equal to the induction 
equation for ideal MHD but with an additional source 
term 



d r PrB = - 



rB 



CT 



(PrB)l 
{(3rB) 2 



(1) 



The index denotes quantities at some initial radius r 
where the dissipation starts. The constant p stands for the 
ratio between the magnetic field component which cannot 
dissipate (described by the ideal MHD induction equation) 
and the total field strength at r = ro. Speaking in terms 
of Poynting flux this means that p 2 is (almost) equal to 
the Poynting flux fraction which does not dissipate, p = 
corresponds to a complete decay of the magnetic field 
while p = 1 means no dissipation at all. 

We derived the functional form of the dissipation time 
scale r in Paper II: 



2nT 2 



1 + u e 



where 



ua 



B r 



B 



V4 



7TW 



V47rr 2 



w 



(2) 



(3) 



is the Alfven 4-vclocity in the comoving frame and w is 
the proper enthalpy density. In regions where the magnetic 
energy density dominates over the matter energy density 
the Alfven velocity is near the speed of light so that the 
square root in (2) is close to 1. This approximation was 
used in Paper II to obtain analytical results. In the present 
numerical study this approximation is not made. 

2.2. Poynting flux and baryon loading 

A very important parameter of our model is the ratio be- 
tween Poynting flux and kinetic energy flux at the initial 
radius r which we denote by 



£ P f,o _ Po^B 2 



CT = 



(4) 



where M is the mass flux per sterad. The outflows of in- 
terest for us are Poynting flux dominated so that <r ^> 1. 
The initial Poynting flux ratio controls not only the initial 
velocity but also the final velocity of the flow as explained 
below. 

How the flow is accelerated from very low velocities 
near the source is not possible to calculate with the ap- 
proach presented since the azimuthal velocity and the 



4 



G. Drenkhahn and H. C. Spruit: Efficient acceieration and radiation in Poynting flux powered GRB outflows 



radial field components cannot be neglected there. We 
assume that magnetocentrifugal acceleration mechanisms 
work there accelerating the flow up to the fast magne- 
tosonic speed as it is the case in stellar winds. For rela- 
tivistic Poynting flux dominated winds we know that the 
typical length scale for this acceleration is on the order 
of the Alfven radius which is in size similar to the light 
radius. If the amount of initially injected thermal energy 
is small (not much larger than the rest mass energy) it is 
converted quickly into kinetic energy so that we can treat 
the flow to be cold again at a few light radii. In the cold 
limit the fast magnetosonic speed equals the Alfven speed. 
The Alfven 4- velocity is a function of the initial Poynting 
flux ratio a only at r — r : ua,o = V^o (Paper II). We 
now take this value as initial 4-velocity for our numerical 
calculations uq — ^/oq which start at ro ~ a few x c/Ct. 

In other GRB studies the baryon loading (mass flux) 
M and the total energy flux L determine the Lorentz fac- 
tor by T = Lj (Mc 2 ). It is often assumed that all the avail- 
able energy is converted into kinetic energy at first. Under 
the same assumption we showed in Paper II that the fi- 
nal Lorentz factor of a dissipating Poynting flux outflow 
is Too = \A + °o ~ vf/ 2 - Hence, the baryon loading M is 
determined from the total energy flux L and the Poynting 
flux ratio cr by M = (7q 3 ^ 2 L/c 2 . It is a matter of taste 
weather one describes an outflow by L, M or L, cto but we 
use the latter to keep the notation of Paper II. 

2.3. The role of electric fields in the flow 



so that one can substitute a c by in the stationary 

(d t B = 0) induction Eq. (5): 



curl (v x B) = — 



\SE\ 



4tt \j\ 



-V 2 B 



(7) 



Using the stationary form of Ampere's law j = c ■ 
curlJB/(47r) to eliminate j and solving (7) for \SE\ gives 
\SE\ as function of v and B only: 



\SE\ = 



| cur IB | |curl(u x B)\ 

c\V 2 B\ 



(8) 



In spherical coordinates and for a radial flow where dIB 
this reads 



5E 



(d r rB)(d r rf3B) 



rd 2 r (rB) 



As fraction of the ideal MHD electric field -E m hd 
that is 



SE 



d r \n\r(3B\ 



d r In | d r rB | 



(9) 

: (3B 

(10) 



This expression is a function of r, (3, B and can be calcu- 
lated numerically. SE/(j3B) takes its maximum at radii 
where the dissipation ceases where it can be of the order 
0.1 depending on the chosen input parameters. For the 
largest part of the flow 5E/((3B) <C 1 so that the use of 
the ideal MHD equation for the evolution of the field is 
justified. The effect of small scale reconnection processes 
is instead taken into account by the decay term in (1). 



In our present study we use the dynamic equations for 
ideal MHD flows and have to make sure that the ideal 
MHD approximation is applicable. In ideal MHD, the elec- 
tric field vanishes in the frame moving with the fluid. But a 
non-vanishing comoving electric field must exist near the 
reconnection centres for field annihilation to take place. 
We assume that the spacial regions occupied by a non- 
vanishing electric field are small and only restricted to the 
reconnection centres. Then we can neglect its influence on 
the dynamic equations on larger scales (see Sect. 2.4). But 
still, the field dissipation produces an extra magnetic field 
gradient also on large scales resulting in X7 2 B ^ in the 
lab frame. We show in this section how the large scale 
electric field can be estimated and that it deviates only 
by a small component SE from the electric field of ideal 
MHD. This component can be neglected in the numerical 
calculations. 

The following calculations are done in the lab frame 
for quantities which vary only on large length scales. We 
start with the induction equation with non- vanishing con- 
ductivity <j c : 



d t B = curl (v x B) + —\7 2 B 

Ohm's law reads 



(5) 



2.4. Dynamic equations 

The conservation equations for mass, energy, momentum 
together with Eq. (1) describing the evolution of the mag- 
netic field determine the proper mass density p, the proper 
internal energy density (excluding the rest mass energy 
density) e, the radial 4-velocity u and the magnetic field 
strength B as functions of radius. In our model the mass, 
energy, and momentum equations read 



d r r pu = 



drT 2 



wTu 



(3B 2 

47T 



P- 



1 



2rp 



(11) 



(12) 



(13) 



(cf. Konigl & Granot 2002; Lyutikov 2001). The variable w 
denotes the proper enthalpy density w = pc 2 + e+p where 
p is the thermal pressure. The thermodynamic quantities 
p,e,p,w are defined in the comoving frame. We assume 
the gas (fully ionised hydrogen) to be ideal with negligible 
heat conduction and an equation of state p — (7— l)e (and 
w = 76 ) where 7 is the adiabatic index. 

The continuity and energy Eqs. (11), (12) are inte- 
grated to give the total mass loss per time and per sterad 



j = a c (E + (3 x B) 



SE 



(6) M 



r 2 upc 



(14) 
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and the total luminosity per sterad 
{rBf 



w 



L = — -TMc + (3c- 
pc z 



47T 



(15) 



where one identifies the kinetic energy flux per sterad 
-kmat = w / (pc 2 )T M c 2 and the Poynting luminosity per 
sterad L pf = (3c(rB) 2 /(4tt). 

2.5. Below and beyond the photosphere 

As long as the medium is optically thick, matter and ra- 
diation can be considered to be in thermal equilibrium. In 
this case, e includes both the thermal particle energy and 
the radiation field energy. The pressure is dominated by 
the radiation so that the adiabatic index is 7 = 4/3. At 
the photosphere radius r p h, where the outflowing material 
becomes optically thin, the radiation decouples from the 
matter and escapes as black body radiation. The pressure 
and internal energy at radii r > r p h is only provided by 
the matter. 

The transition from optically thick to optically thin 
conditions is sharp, in practical flow models, and we sim- 
plify the computations here by treating it as a disconti- 
nuity. Its location is in principle found by integrating the 
optical depth into the flow. As in stellar atmospheres and 
winds, a fair approximation for its location is the point 
where the mean free path of a photon equals the density 
scale height. In our case, this is of the order of the distance 
r from the source. At the photospheric temperatures en- 
countered (a few keV) the dominant opacity is electron 
scattering. Taking into account the Lorentz transforma- 
tion of the mean free path to the rest frame, the photo- 
spheric radius is thus given by r p h = T/(nThp) where KTh 
is the Thomson scattering opacity and p the density in the 
comoving frame. 

For r > r p h the thermal pressure is only supplied by 
the plasma and the adiabatic index 7 depends on the tem- 
perature. For non-relativistic temperature kT <§C m c c 2 we 
have 7 = 5/3 but for hotter medium the electrons become 
relativistic which lowers 7. Since fast radiative cooling is 
a model assumption the matter in the optically thin part 
stays cold enough so that 7 = 5/3 is valid there. The 
validity is checked during the numerical computations. 

2.6. Radiative loss 

In the optically thin regime energy and momentum from 
the dissipating magnetic field is transfered into radiative 
form. The radiation escapes and does not interact with 
the matter. Let A be the emissivity of the medium in 
the comoving frame, that is the energy which is radiated 
away per unit time and per unit volume. If the emission is 
isotropic in the comoving frame the energy and momen- 
tum Eqs. (12), (13) including the radiative loss terms are 



d r r 2 ( wTu + 

\ 47T 



9 A 

-r 2 r- 

c 



(16) 



d r r 



B 2 \ A 
2 1 wu 2 +p + (l + f3 2 ) — ) =2rp- r 2 u— 

' 87T / c 



(17) 



(Granot & Konigl 2001). 

The importance of the cooling term depends on the 
cooling time scale. If it is short, the matter stays cold (gas 
pressure negligible) during the dissipation process. In this 
limit, all the dissipated energy is locally radiated away. 
Synchrotron emission is a plausible fast cooling process. 
It is particularly effective in our model, because the mag- 
netic field strengths are high in a Poynting flux dominated 
outflow. With (30) one derives a typical field strength 



B <7 ■ 10 7 G 



7-1/2-1 

^50 '13 



In the comoving frame this is 



B, 



co, 3 



£7 



in 2 . j 1 ! 2 ^^- 

iU ^50 r 13 1 2 



(18) 



(19) 



The distance travelled by the medium in one cooling time 
in the lab frame is r coo i = ct syD T. The synchrotron cooling 
time scale in the comoving frame is i syn = 6 s • T~2B~ 2 3 
(Daigne & Mochkovitch 1998) where T e = 100r Cj2 is the 
Lorentz factor of the radiating fast electrons and B co = 
1000 G • -Bco, 3 is the comoving magnetic field strength. In 
units of the expansion length scale of the flow, r, the cool- 
ing time is 



^cool 



4-10" 



• risL^TlT-l « 1 



(20) 



This shows that synchrotron cooling is fast for these fidu- 
cial parameters. Though the simplifying assumption of 
fast cooling in the optically thin regime is thus justified, 
we have kept an ad hoc cooling term in the calculations 
to ease the numerical treatment. 

The form of this cooling term used is 



A = k 



ecu 



(21) 



where k is an adjustable cooling length parameter. The 
cooling length is the distance by which the matter travels 
outward while the internal energy e is lost. We have used 
k = 10 4 so that the cooling length is the distance 10~ 4 r. 
This is only a small fraction of the expansion length scale 
r and thus qualifies for the description of a fast cooling 
flow. 

Because of the fast cooling the temperature is always 
very low (kT <C m e c 2 ) and the equation of state of the 
gas is that of a non-relativistic fully ionised gas, 7 = 5/3 
for r > r ph . 



2.7. Computational method 

We choose 



/ r 2 pc 2 \ 
r 2 e 
u 

\ ^M= I 



Q = 



(22) 
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to be the vector of primitive variables. There are no prin- 
ciple reasons against taking e.g. (p, e, u, B) instead but the 
use of (22) simplifies the following analytical expressions 
a bit. The set of Eqs. (11), (16), (17), (1) can be written 
in matrix form 



A ■ d r q = s 
with the matrix 

( u 

uT juT 

u 2 7 r 2 -i 

\ 



(23) 



q(ro; (Jo, L). The value of tq is unimportant for the flow at 
larger radii (Paper II) and we set r = 3 • 10 7 cm for all our 
calculations without introducing an additional restriction. 
In total the model's input parameter space is effectively 
made up of cr , L, yU, eO. 

Eq. (23) is a ordinary differential equation and can be 
solved numerically with common software packages. The 
integration proceeds stepwise from r = r until the photo- 
sphere is reached, where the mean free path for the pho- 
tons is equal to the radius r. The photosphere must be 
treated in a special way because it is a discontinuity where 
the radiation decouples from the matter part. 



r 2 pc 2 



r 2 wT (l + (i 2 
2r 2 wu 



f p-3 

r 3 4n 

rB p-3 



rB 

/47T 





2(3 



(24) 



and the source term vector 

( A 

_ r 2p7A 

2r(7 - l)e - r 2 u^ 



V 



1 rB 

CT JTlX 



1 ,.2 (»-£Qo 



(25) 



The elements of A and s are functions of r, q and the 
constant model parameters ao, H, efl, Vq. For the cooling 
term A we use the expression 



^ _ i for r < r ph 

10 4 ecu/ r for r > r p h 



(26) 



the dissipation time scale (2), and the adiabatic index 

(27) 



J 4/3 for r < r ph 
^ \ 5/3 for r > r ph 



2.7.1. Boundary conditions and solution process 

To initialise the solver at the initial radius ro we need to 
determine the vector of primitive variables qo = q(ro). 
The flow starts with Alfven velocity go, 3 = u o = and 
is cold q .2 = r%e = 0. By solving (4), (14), (15) for M 
we obtain 



r 2 puc = M 



and thus 



c 2 (a + 1) 3/2 



9o,i = fopoc 2 = 



L/c 



(28) 



(29) 



These three equations can also be solved for ro-Bo to give 
finally 



90,4 



47T 



CO 



(70 + 1 



1/4 



(30) 



The initial vector depends only on the initial radius, the 
initial Poynting flux ratio and the total luminosity: qo = 



2.7.2. Transition at the photosphere of the flow 

At the photosphere, the equation of state changes from 
one dominated by radiation to one dominated by the gas 
pressure. To connect the two, the radiation emitted at the 
photosphere has to be taken into account. The amount of 
energy involved can be substantial, and appears as an (ap- 
proximate) black body component in the GRB spectrum. 
It depends on the temperature of the photosphere. 

The temperature at the photosphere is kT <C m e c 2 for 
all used parameter values so that pairs can be neglected. 
The photosphere is then simply determined by 



r 



rpKTh 



1 



(31) 



At the photosphere one has to subtract the energy and 
momentum which is carried away by the decoupled radia- 
tion. To calculate these quantities one needs the temper- 
ature at the photosphere. 

The dimensionless temperature 6 = kT / (m e c 2 ) in the 
optically thick region is given by the solution of 



m c 2 7T 2 m c c 2 4 



(32) 



where A is the electron Compton wave length. 

At the photosphere we calculate the temperature # p h 
and subtract the radiation energy density of a black body 



ebb = 



7T 2 m 



15 



r 2 

^3~ V 



(33) 



from the total energy density: e = e — ebb- The integration 
proceeds with an adiabatic index of 7 = 5/3. The temper- 
ature # p h is the temperature of the emitted black-body 
radiation which has a luminosity per sterad of 



' 



for r < r ph 



r ph| e bbW p hr ph c for r > r ph 



(34) 



The integration continuous until the dissipation ceases. 
There, the luminosity of emitted non-thermal radiation is 
determined by 



L n t — L — L p f — L mat — L 



bb 



(35) 
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2.7.3. Modifications to the black-body radiation 

We have assumed in the above that the radiation emitted 
at the photosphere is perfect black body radiation which 
means that the photons are created at the photosphere. 
But if scattering dominates (scattering coefficient is of the 
same order or greater than absorption coefficient) the pho- 
tons from the photosphere are produced at smaller radii 
(and at different temperatures) and undergo many scat- 
terings until they escape to infinity. 

The number density of the photons is determined 
at the radius where they are created. The relevant cre- 
ation processes are free-free and Synchrotron emission. 
Depending on the degree of Comptonisation the emergent 
photon spectrum will be a modified black body or a Wien 
spectrum (Rybicki & Lightman 1979). The typical photon 
energy is then given to the photons by Comptonisation at 
some radius larger than the creation radius but smaller 
than the photospheric radius. We thus have three char- 
acteristic radii for the radiation process: the photon cre- 
ation radius where the number density is determined, the 
Comptonisation radius determining the typical photon en- 
ergy and the photosphere radius where the radiation de- 
couples from the matter. 

We can speculate about the spectrum of the emer- 
gent radiation if scattering processes are included. In ther- 
mal equilibrium the radiative pressure-temperature rela- 
tion reads p a T 4 . But if the photon number does not 
depend on temperature any more the relation changes to 
p oc T leading to an increase in temperature (similar to 
the one seen at the photosphere in Fig lc). The photon 
production rate also rises until some equilibrium tempera- 
ture is reached. Thus radiation and temperature cannot be 
determined independently. Compared to the black body 
temperatures at the photosphere as used in the model 
presented the real temperatures may thus be somewhat 
higher. The number of photons might be less since they 
originate from a smaller area. Because the absorption and 
emission processes are highly frequency dependent pho- 
tons of different energy might not be produced at the 
same radius and with the same bulk Lorentz factor. This 
could lead to a broader spectrum though reprocessing by 
Comptonisation may also play a role. 

A detailed modelling of the emergent radiation at the 
photosphere requires a consistent but rather complicated 
extension of the model presented here. We have to assume 
that the simplified black body treatment at least give rea- 
sonable estimates for the total energy though the spectral 
shape might differ in reality. The task to determine the 
spectrum must be postponed to a more detailed investi- 
gation in the future. 

2.7.4. Another transition radius 

The rate of dissipation of the magnetic free energy starts 
out fast, but as the field strength decreases it become 
slower than the expansion time scale. There is thus a char- 
acteristic radius, which we call here the saturation radius 



r sr , beyond which magnetic dissipation effectively stops. 
Since the acceleration of the flow is intimately connected 
with the dissipation, this is also the radius where the flow 
reaches its terminal speed. In Paper II we have shown that 
r sr is of the order 



(36) 



while the terminal Lorentz factor Too is of the order 

T x = (1 - p 2 ) 4 /2 . (37) 

A simple approximation for the dependence of the Lorentz 
factor on distance then turns out to be 



Loc (r/r sr ) 1/3 for r < r s 
Loo for r > r s 



(38) 



2.7.5. Solution examples 

As an example, Fig. 1 shows the result of a numerical in- 
tegration. The result can be compared with the analytical 
approximation derived in Paper II. The analytic estimate 
gives a fair representation of the full results, though it 
overestimates the terminal Lorentz factor somewhat. 

The luminosity carried by the medium L mat — 
w/(pc 2 )TMc 2 is made up of the kinetic Lkin = TMc 2 and 
the thermal part L th = (w/(pc 2 ) - l)rMc 2 . The fact that 
the fluid is dominated by the pressure and energy density 
of the radiation in the optically thick region can be seen 
at the photosphere. Fig. lc displays that a major part of 
the thermal energy flux L t h is made up from the radia- 
tion component which is released as black body radiation 
at the photosphere. This explains why the thermal energy 
flux L t h in Fig. lb nearly coincides with the black body 
luminosity Lbb- 

Outside the photosphere T (and therefore £kin) be- 
comes smaller than the analytical estimate. The non- 
thermal radiation flux component L Dt rises quickly and 
the dissipated energy is efficiently converted into radia- 
tion. The fractions of the total luminosity converted to 
kinetic, thermal, and non-thermal energy in this example 
are 54%, 5%, and 41%, respectively. This demonstrates 
the high efficiency of the magnetic dissipation process in 
accelerating the flow and the production of non-thermal 
radiation. 

In Paper II we neglected the thermal pressure and one 
could argue that the acceleration obtained in the present 
study should be somewhat greater due to the additional 
effects of the finite thermal pressure gradient of the mat- 
ter and radiation component. To compare the magnetic 
and thermal contribution to the acceleration we regard 
the momentum eq. (13) which reads in non-conservative 
form 



^rd r r 2 wu 2 + d r p=- (j x B) 



d r (rB) 2 
8irr 2 



(39) 



The ratio between the thermal pressure gradient d r p and 
the magnetic force density [9 r (ri3) 2 ] / (87rr 2 ) is plotted 
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Fig. 1. Solution for L = 10 50 ergs -1 sterad~\ a = 100 (cor- 
responding to an initial Lorentz factor of 10), ef2 = 10 3 s _1 , 
H = 0. The vertical dotted line indicates the photosphere ra- 
dius, (a) Lorentz factor where the dotted line represents the 
analytical approximation derived in Paper II. (b) Various en- 
ergy fluxes per sterad as labelled in the legend. Indices de- 
note the Poynting flux (pf), fcmetic, i/iermal, black-body (bb) 
and non-thermal (nt) components, (c) Normalised temperature 
= kT/(m c c 2 ). The matter is rather cold and pairs do not play 
a role since 6 <C 1. The temperature jump at the photosphere 
is no real discontinuity but rather a rapid but steady change. 
This is the result of the abrupt change in the equation of state 
from a radiation dominated to a matter dominated fluid, (d) 
Ratio between the thermal pressure gradient and the magnetic 
force density in the flow. The acceleration is completely deter- 
mined by the magnetic field at all radii. 



3. Shortest time scales 

Though the model presented here is stationary and thus 
does not describe the variability of GRBs, we can check 
whether the physical model on which it is based is compat- 
ible with the observed millisecond variations. If the out- 
flow contains inhomogeneities, say, regions where the re- 
connection proceeds faster or unsteady the flow will show 
more ore less bright patches. In this section we show that 
the emission time of these patches will be short enough 
for GRBs. 

The time it takes for the magnetic field to dissipate 
sets the limit for the shortest emission duration of a sin- 
gle patch in the flow. The time interval dt — dr/c in which 
a flow element moves outward by dr is Doppler boosted to 
the observed time interval di bs = (1 — /3) dt w dr/(2cT 2 ). 
Let us assume that the more inhomogeneous flow can 
still be described approximately by the stationary solu- 
tion. Using (38) and integrating the observed time from 
r = . . . r sr gives 

tabs — t^-^ = 1-6 • KT 3 s • (e_ 1 4 )- 1 . (40) 

If the flow consists of shells with differently strong or 
fast reconnection/emission it can produce observable vari- 
ations on time scales of the order of milliseconds. 

To see if small patches can contribute a significant vari- 
ability, we note that their lateral size (perpendicular to 
the flow) is likely to be of the same order as A = 2-kc/SI. 
This is the typical size of the small scale field inhomo- 
geneities and the reconnection centres. Due to the short 
cooling time the emitted radiation originates from close to 
the reconnection centres so that A is also the lateral size of 
the bright patches. The overall expansion in the flow does 
not change the size of the patches if the expansion speed 
in the comoving frame is always small compared to the 
Alfven speed w c which is the typical speed with which 
the magnetic field can reorganise itself. This is the case for 
A < r/r being always true at the radii of interest. Thus 
these patches will stay of the same size A =const. and do 
not scale with radius. A is quite small compared with the 
radius of the photosphere. The solid angle from which the 
observed radiation is emitted is of the order 1/r 2 . Hence 
there arc n = (r/XT) 2 = 3 • 10 11 • 51|r^" 2 r 2 3 reconnection 
patches contributing at any moment to the observed ra- 
diation, and the maximum variability amplitude expected 
is thus y/n/n w 2 • 10~ 6 . The observed time scales are 
not produced locally in the flow but must be due to a 
variability of the central engine. 



in Fig. Id for our chosen example parameter set. At the 
photosphere the ratio decreases since the radiation compo- 
nent drops out. But even if radiation is included the total 
thermal pressure gradient is not important compared to 
the magnetic force density which determines the dynamics 
of the flow. Thus the radiation pressure can be neglected 
above the photosphere without influencing the dynamics 
of the flow. 



4. Parameter study 

The parameters of the model are the initial Poynting flux 
to kinetic energy flux ratio a , the total luminosity per 
sterad L, the fraction of non-dissipatable magnetic field 
/i, and the measure for the reconnection rate ef2. In this 
section we explore the dependence of the solutions on these 
parameters, by plotting values of the physical quantities 
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Fig. 2. Some physical quantities at the photosphere of the flow 
as function of the initial Poynting flux ratio ao for different 
values of the total luminosity per sterad. The values efi = 
1CF 3 s _1 , fx = where used, (a) Temperature in the comoving 
frame, (b) Lorentz factor at the photosphere, (c) Radius of the 
photosphere. Straight lines in panel (b) + (c) at 200 < ao < 
1000 indicate the analytical solutions from (41), (42). Thick 
dots correspond to ao = <7o,br defined in (43). 



at the photosphere and their asymptotic values at large 
distances. 



4.1. The photosphere 

Figure 2 shows the temperature, Lorentz factor and radius 
of the flow at the photosphere as functions of the initial 
Poynting flux ration ao- In all plots of the three quanti- 
ties one can identify a break at a certain value of ao- This 
break can be understood in terms of the saturation radius 
r sr mentioned in Sect. 2.7.5. The reconnection yields only 
little energy beyond r sr since the largest part of the free 
magnetic energy is already used up before. The location 
of this radius relative to the photosphere determines the 
basic properties of the results. In order to interpret the 
numerical results, we recall here some the results of the 
analytic model derived in Paper II. In this model the ra- 



dius of the photosphere r p h and the Lorentz factor at the 
photosphere are given by 



1.05 • 10 n cm 



€-1^4 



0.5 



-2/5 



r 3/5 -3/2 
-^50 °0.2 ' 



Tph ~ M ph = 119 • 



1-^2 
£-1^4 [ - - ) L 50 



1/5 



(41) 
(42) 



By equating r p h and r sr of (41), (36), we find the value of 
ao where the break in the parameter dependences occurs: 



co,br — 39 • (e^ifl^L^oY 



2/15 ( 1 - M 2 
0.5 



-8/15 



(43) 



At this value of the baryon loading, the dissipation of mag- 
netic energy ceases around the photosphere. For ao < ao.hr 
dissipation occurs mostly inside the photosphere, and the 
radiation is dominated by a black body component. 

We note that the analytic expression (41) for r p h has 
been derived under the assumption that most of the dissi- 
pation occurs outside the photosphere. It is still accurate 
enough, however, for the estimate (43) that we use to in- 
terpret the numerical results. The asymptotic validity of 
the analytic value (42) of the terminal Lorentz factor 
for large ao is shown in Fig. 2b. 

4.2. Limits on thermal and non-thermal radiation 

The model yields the luminosity per sterad for both the 
black-body radiation from the photosphere and the non- 
thermal radiation. In this section we investigate how these 
radiation components behave as function of the model pa- 
rameters and what observed temperatures are expected for 
the black-body component. 

Figure 3b displays the luminosities of the thermal and 
non-thermal radiation components Lbb and L nt as frac- 
tion of the total luminosity L. At very low co-values both 
components are very small. In that case the energy is re- 
leased far below the photosphere and is converted into ki- 
netic energy. The same happens in 'dirty fireball' models 
where the central engine injects thermal energy into the 
matter near the source (Shemi & Piran 1990; Paczynski 
1990). This also leads to an almost complete conversion 
into kinetic energy. 

The black-body radiation shows its maximum if the 
dissipation ceases right at the photosphere so that r p h = 
r sr . This corresponds to ao = a .bi from (43). This analyt- 
ical estimate for er ,br might be not coincide exactly with 
the maximum of the numerically obtained Lbb / L-curves in 
Fig. 3b since some simplifications were used in the deriva- 
tion of (43). Though, we take Lbb, max = Lbb(co,br) to be 
the maximal black-body luminosity to simplify the treat- 
ment in the following. 

From the graphs in Fig. 3b one might guess that the 
maximal value of L^b/L docs not depend on L. Indeed, 
it turns out that the maximal fraction of the black-body 



10 



G. Drenkhahn and H. C. Spruit: Efficient acceieration and radiation in Poynting flux powered GRB outflows 




Fig. 3. Terminai Lorentz factor and radiation flux ratios as 
function of the initial Poynting flux ratio <to f° r different values 
of the total luminosity per sterad. The values efi = 10 _3 s _1 
and p — where used, (a) Terminal Lorentz factor Too. 
The dotted line correspond to the analytical estimates from 
Paper If where ao 3> 1 and no radiative losses were assumed, 
(b) Ratio between black-body and total luminosity Lbb /L and 
ratio between non-thermal and total luminosity L nt /L. At the 
location of the thick dots, the model parameters are such that 
the magnetic dissipation ceases to be effective near the photo- 
sphere [cf. (43)]. 
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Fig. 4. Maximal thermal and non-thermal emission as function 
of efi (a measure for the dissipation rate), (a) Ratio of maxi- 
mal thermal emission to dissipatable luminosity Lbb,max/[(1 — 
li 2 )L]. The used values for Lbb, max correspond to the thick 
dots in Fig. 3b. (b) Ratio of maximal non-thermal emission 
to dissipatable luminosity L n t,max/[(1 — /i 2 )L] determined at 
a = 10 3 . 



luminosity to the dissipatable luminosity is almost a con- 
stant. In the Poynting flux dominated wind the initial 
Poynting flux luminosity is almost equal to the total lumi- 
nosity ipf.o ~ L. The fraction of Poynting flux which can- 
not dissipate by reconnection was parameterised by p? so 
that dissipatable luminosity is (1 — /i 2 )L p f.o ~ (1 — H 2 )L. If 
we plot Lbb. max/ [(1 ^ ^ 2 )L] as a function of the two other 
model parameters p, efi we see that its value is around 
0.17 ± 0.03 as displayed in Figure 4a. Thus, the energy in 
black-body radiation is always less than 20% of the total 
releasable magnetic energy. 

Figure 3b shows that the fraction of the total lumi- 
nosity emitted as non-thermal radiation has a maximum 
value for large ao of about 50%, independent of the lumi- 
nosity itself. In this limit almost all the dissipation takes 
place outside the photosphere and L^/L is negligible. 

Figure 4b shows that L nt;lnax /[(1 — /x 2 )L] is always very 
close to 0.5. The maximal radiation efficiency occurs in the 
extreme Poynting flux dominated limit. A fast radiation 
mechanism converts half of the free magnetic energy into 
non-thermal radiation. 

The dissipatable energy flux (1 — /i 2 )L, not the total 
Poynting flux L p f m L in general, is the energy reser- 
voir from which the radiation and kinetic energy is fed. 
Nevertheless, one needs to know the total Poynting flux in 
order to determine the absolute magnetic field strength in 
the medium to investigate the physical emission process. 
Since this is not needed here we could restrict our study 
to the case \x = in the largest part of this paper. For 
any combination of /z, L one finds p! = 0, V = (1 — p 2 )L 
which yield an outflow of equal dissipatable Poynting flux 
and thus equal emission and kinetic energy. Therefore, the 
setting [i — has not introduced an additional limitation 
in the generality of our investigation. 

4.3. Limits on the terminal Lorentz factor 

In Paper II we derived the terminal Lorentz factor for a 
complete conversion of Poynting flux into kinetic energy 
flux. For a Poynting flux dominated flow this terminal 
Lorentz factor is T x = (1 — /i 2 )^ 3 / 2 . Since we found that 
the radiative losses can be as large as 50% of the total 
luminosity we can now determine upper and lower limits 
for the terminal Lorentz factor: 

± (l - M >S /a < Too < (l - M 2 )^ • (44) 

4.4. Possible variability 

The stationary treatment in our study docs not yield 
any variability by definition. But we can speculate about 
the outcome of quasi-stationary changes in one or more 
model parameters. The non-thermal luminosity L nt de- 
pends most strongest on ao at moderate values. This is 
seen in Fig. 3b where the L nt /L graphs rise quickly to 
the maximal values in an values around 20 ^ <tq 100. 
A variation of ao in this moderately large interval might 
produce some kind of on-off behaviour of the non-thermal 
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Fig. 5. Thermal to non-thermal flux ratio Lbb/int and red- 
shift corrected temperature of thermal radiation for various 
model parameters. Calculations for each panel are done for 
fixed values of fj, = and efi as indicated. Solid lines corre- 
spond to parameter sets with equal total luminosity L while 
dotted lines correspond to equal initial Poynting flux ratio <7 () . 



luminosity and the large variability in GRB light curves. 
Extending the model to include time dependence remains 
an interesting investigation for the future. 

4.5. Observable quantities 

The results show that a thermal component in the emis- 
sion is expected at higher baryon loading values. In a lim- 
ited range of baryon loading (10 ^ a < 70), the model 
predicts a GRB with a significant thermal component. 
This property can be used observationally as a test of the 
model, or as a diagnostic of the GRB outflow. 

Our model GRBs can be represented in a diagram 
of showing the ratio Lbb/^nt as a function of the black 
body temperature. This is shown in Fig. 5. The photo- 
spheric (black body) temperature T b s shown is the value 
as observed in the rest frame of the GRB host. It is re- 
lated to the temperature in a comoving frame Tb b by 
T bs = Tphlbb- 



The single panels of Fig. 5 show lines for which 3 of the 
4 model parameters are fixed while L or <r is varied. An 
increase in the initial Poynting flux ratio o$ increases the 
temperature but decreases the black-body component due 
to a smaller photosphere radius (see Fig. 2c). An increase 
of the total luminosity L results in a larger photosphere 
temperature and a larger black-body component though 
the dependence of Lbb / L nt and fcT Q b s on L is much weaker 
than on <jq. 

From Fig. 5 one finds that only low Poynting flux ratios 
cr ^ 200 lead to a significant fraction of thermal radiation 
Lbb/L n t <; 0.1. The predicted black body temperatures 
can range all the way from about 5 to 100 kcV. These are 
the temperatures as observed in the frame of the GRB 
host, so a redshift has to be known in order to use these 
predictions diagnostically. 

Thermal components have indeed been observed in 
GRB spectra. Preece (2002) reports on a thermal compo- 
nent in the spectra of GRB 970111 within the first « 10 s 
after the trigger. In this time interval the temperature of 
the black body component was observed to vary between 
45 and 75 keV, and the ratio of thermal to non-thermal 
flux was of the order unity. After this initial phase the non- 
thermal component started to dominate. In terms of our 
model, these observations indicate that this GRB started 
with a moderate baryon loading, which then decreased in 
the course of the burst. Since a redshift of this burst has 
not been determined, a more detailed comparison with the 
model can not be made. 

4.6. Connection with X-ray flashes 

X-ray flashes are fast X-ray transients which are not de- 
tected in the 7-ray band 40-700 kcV of BcppoSAX (Heise 
ct al. 2001; Heise & in 't Zand 2002). In our model this 
spectral characteristic can be explained by an outflow of 
low <7o where the thermal radiation dominates over the 
non-thermal component. The relevant region in Fig. 5 
would be at Lbb/L nt <; 1. The temperatures predicted 
by the model for this range of black body luminosity is 
^ 30 keV, which is quite compatible with the observations. 

5. Summary and discussion 

A magnetised and rotating central engine of a GRB pro- 
duces a Poynting flux outflow. Any non-axisymmctry of 
the magnetic field leads to small scale (wave-like) varia- 
tions in the electromagnetic field carrying energy outward. 
We have assumed there that these small scale irregular- 
ities are subject to rapid reconnection, governed by the 
Alfven speed, as observed in other astrophysical settings 
and in numerical simulations. Thus, the magnetic field can 
rearrange itself to a energetically favourable configuration 
and releases its free energy stored in the small scale field 
variations. 

The release of free magnetic energy proceeds with a 
rate determined by the length scale of the field variation 
and the local Alfven speed of the plasma. The magnetic 
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field acts as energy reservoir carried with the matter which 
transfers its energy to the matter continously. The decay 
of the magnetic field at the same time causes an outward 
gradient of the magnetic pressure. This causes a signifi- 
cant part of the Poynting flux to be converted into kinetic 
energy The other part of the free energy is converted into 
heat. In the optically thick region of the flow a thermal en- 
ergy gradient promotes adiabatic expansion and the con- 
version of thermal energy to kinetic energy Thus, at small 
radii in the optically thick region, almost all of the dissi- 
pated magnetic energy gets converted into kinetic energy. 

When the flow becomes optically thin at the photo- 
sphere the thermal radiation energy escapes to infinity. 
This radiation resembles a black body spectrum if the pho- 
tons are created near the photosphere. But if scattering 
dominates the spectral shape will differ and might look like 
e.g. a modified black body or a Wien spectrum (Rybicki 
& Lightman 1979). A much more detailed treatment is 
necessary to determine the radius at which the radiation 
is produced and what spectral shape it has. At this stage 
we are only interested in the total energetics and the char- 
acteristic radiation temperature and we assume that the 
simplified black body treatment give sufficiently precise 
estimates. 

To compute the radiation spectrum from the optically 
thin region more detailed physics is needed which is be- 
yond the scope of this paper. We have instead assumed 
that the reconnection process under optically thin condi- 
tions maintains a significant population of energetic elec- 
trons, which then radiate synchrotron radiation in much 
the same way as in the standard internal shock model. An 
advantage of our model is that the magnetic field needed 
for the synchrotron radiation is a natural part of the flow 
model itself (see also Paper I). The central difference of 
the model presented to standard internal shock models for 
GRBs is that radiation stems from the local dissipation of 
magnetic energy and not from shock conversion of kinetic 
energy. Therefore, one does not need an extremely variable 
central engine to obtain an acceptable radiation efficiency 
(Beloborodov 2000; Kobayashi & Sari 2001). 

The continuous character of the energy release leads 
to a slower acceleration of the flow compared to the clas- 
sical fireball scenario. In a fireball the energy is injected 
abruptly as thermal energy. This leads to an rapid accel- 
eration where the Lorentz factor is linear to the source 
distance r cx r. In our model the release of magnetic en- 
ergy leads to T cx r 1 / 3 in the optically thick region. 

An important model parameter is the ratio between 
Poynting flux and kinetic energy flux <r at some initial 
radius r . This parameter controls the baryon loading 
parameter in a sense that high values correspond to a 
low baryon loading. The value of a decides how much 
of the Poynting flux energy gets converted into kinetic 
energy, black-body radiation and non-thermal radiation. 
The three other parameters of the model are the total lu- 
minosity per sterad, the fraction of dissipatable Poynting 
flux and the reconnection speed. There are three intervals 
of <Jo values in which the characteristics of the flow is sig- 



nificantly different. At very low values (a ^ 10) most of 
the energy gets converted into kinetic energy The mag- 
netic energy gets released in the optically thick part. This 
case is similar to dirty fireball models in which the cen- 
tral engine injects thermal energy into the matter near the 
central engine (Shemi & Piran 1990; Paczyhski 1990). The 
matter is already cold when it reaches the photosphere at 
large radii and there is no more free magnetic energy avail- 
able to power the non-thermal radiation. The burst energy 
can only power an afterglow by an external shock. 

An intermediate Poynting flux ratio of Co ~ 100 causes 
the release of a considerable amount of energy near the 
photosphere and the thermal emission is non-negligible. 
The black-body component becomes maximal if the ra- 
dius where the dissipation ceases coincides with the pho- 
tosphere. Then, s=a 17% of the dissipatable magnetic en- 
ergy gets converted into black-body radiation while the 
rest ends up in kinetic energy. 

At very high <r ^ 300 values the radius of the pho- 
tosphere is small and almost all of the dissipation takes 
place in the optically thin region. The dissipated energy 
gets equally distributed among the non-thermal radiation 
and the kinetic luminosity of the flow. At low baryon load- 
ing, and for a purely non-axisymmctric magnetic field, al- 
most exactly 50% of the the Poynting is converted into 
kinetic energy and 50% into non-thermal radiation. If a 
substantial part of the magnetic field is axisymmetric, the 
Poynting flux associated with it does not dissipate, and 
instead is expected to show up as afterglow emission. 

Assuming that regular GRBs have large Co this find- 
ing predicts that the energy of the afterglow (fed by the 
kinetic energy of the flow) is comparable to the energy 
of the the prompt emission. Beaming effects change this 
picture if the outflow consists of a sufficiently narrow jet. 
The energy in the afterglow will be weaker because af- 
ter the flow has decelerated the radiation is spread over 
a larger solid angle compared to the highly beamed ini- 
tial radiation (cf. the light curve break discussion in e.g. 
Ghisellini 2001). Therfore the prompt emission might be 
more luminous than the afterglow luminosity if the jets 
points towards us. 

Besides the luminosity of the black-body radiation the 
model yields the temperature of this radiation. We find 
that the unredshiftcd observable temperature is 5 keV ^ 
fc?obs ^ 100 keV for our fiducial GRB/X-ray burst pa- 
rameters. The model produces a rather constant Lorentz 
factor at the photosphere so that large variations of model 
parameters result in only a small temperature spread. We 
cannot make a clear statement about the contribution of 
the non-thermal component to emission in the quoted en- 
ergy range because we do not know the radiation mech- 
anism. If the thermal component at its maximum is of 
the same order or greater than the non-thermal compo- 
nent a feature should be present in the spectrum. In fact, 
there exist observations of excess emission in the low en- 
ergy range (« 1-5 keV) for some GRBs (Strohmayer et al. 
1998; Preece et al. 1996). Recent investigations by Preece 
(2002) show clearly a strong thermal component in the 
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during the first 10 s of GRB 970111 with temperatures of 
45-75 keV. These numbers are in agreement with some of 
our fiducial GRB model parameter values. Observations 
of this kind will enable us to determine the model param- 
eters and even there time-dependence during the burst. 

For initial ratios between Poynting flux and kinetic 
energy flux oo ~ 40 the black-body radiation component 
dominates over the non-thermal component. The radia- 
tion efficiency is lower than in the Co <^ 300 case because 
only ^ 17% of the total luminosity can be converted into 
radiation. We speculated in Paper II about the possibility 
that these low-cro outflows could be identified with X-ray 
flashes observed by BeppoSAX (Heise et al. 2001; Heise & 
in 't Zand 2002). Because no dissipation takes place out- 
side of the photosphere there is no non-thermal emission 
in the 7-ray range > 40keV. The present study showed 
that the thermal emission has an (unredshifted) observ- 
able temperature of ^ 30 keV which agrees with the ob- 
servations of X-ray flashes. Heise et al. (2001) speculated 
that X-ray flashes are bursts with high mass loading. This 
is also true in our model since high mass loading corre- 
sponds to low do values. 

The hypothesis that the thermal radiation from a low- 
cro outflow produces a X-ray flash may be checked by fu- 
ture observation. The model predicts a lower radiation ef- 
ficiency of ^ 17% so that most of the dissipated energy in 
the outflow goes into kinetic form. It will be converted into 
radiation in the external shock of the afterglow. The after- 
glow of an X-ray flash should be more luminous than the 
prompt emission if jet-effects do not interfere too much. 

A stationary approximation for the flow is used in this 
paper. If the central engine operates intermittently inter- 
nal shocks could occur in the magnetised outflow. One 
could also imagine that not the total luminosity changes 
with time but that the other wind parameters like the 
mass loading and with it <r varies, cto has a strong in- 
fluence on the non-thermal luminosity L nt and the black- 
body luminosity Lbb- A time- varying ctq around interme- 
diate values leads certainly to large modulations in the 
non-thermal light curve. The physical model should be 
extended to include time dependent model parameters to 
investigate their effect on the light curve. 

The minimal observed variability of GRBs is around 
a millisecond. Any process producing the emission must 
thcrforc be fast enough to account for this limit. The sta- 
tionary model predicts that the reconnection in the co- 
moving frame lasts for approximately 1 millisecond. The 
effects for Doppler shift and relativistic time dilation al- 
most cancel so that one observes almost the time in the co- 
moving frame. The millisecond variability is therfore com- 
patible with the reconnection model. If the reconnection in 
the flow is not smoothly distributed but patchy we would 
expect to see a peak for the emission coming from one 
of the patches where reconnection takes place. Our model 
is nevertheless applicable because we only need the over- 
all reconnection rate, the average over small length scales, 
which is responsible for the global flow dynamics. 
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